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ABSTRACT 

We analyze the scenario of baryogenesis through leptogenesis induced by the out-of- 
equihbrium decays of heavy Majorana neutrinos and pay special attention to CP violation. 
Extending a recently proposed resummation formalism for two-fermion mixing to decay 
amplitudes, we calculate the resonant phenomenon of CP violation due to the mixing of 
two nearly degenerate heavy Majorana neutrinos. Solving numerically the relevant Boltz- 
mann equations, we find that the isosinglet Majorana mass may range from 1 TeV up to 
the grand unification scale, depending on the mechanism of CP violation and/or the flavour 
structure of the neutrino mass matrix assumed. Finite temperature effects and possible 
constraints from the electric dipole moment of electron and other low-energy experiments 
are briefly discussed. 
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1 Introduction 



Based on the assumption that the Universe was created initially in a symmetric state with a 
vanishing baryon number B, Sakharov [|l|] derived the three known necessary conditions that 
may explain the small baryon-to-photon ratio of number densities n^/n^ = (4 — 7) x 10^^°, 
which is found by present observations. The first necessary ingredient is the existence of B- 
violating interactions. With the advent of grand unified theories (GUT's), this requirement 
can naturally be fulfilled at very high-energy scales 0, through the decay of super- heavy 
bosons with masses near to the grand unification scale Mx ~ 10^^ GeV. However, such a 
solution to the baryon asymmetry in the Universe (BAU) faces some difficulties. In fact, 
Sakharov's second requirement for generating the BAU prescribes that, by the same token, 
the 5- violating interactions should violate the discrete symmetries of charge conjugation 
(C) and that resulting from the combined action of charge and parity transformations 
(CP). One major drawback of the solution suggested is that minimal scenarios of grand 
unification generally predict very small CP violation, since it occurs at very high orders in 
perturbation theory. Therefore, one has to rely on no-minimal representations of GUT's 
in order to obtain appreciable CP violation. Furthermore, experiments on the stability of 
the proton put tight constraints on the masses of the GUT bosons mediating B violation 
and their couplings to the matter. 

The most severe limits on scenarios for baryogenesis, however, come from Sakharov's 
last requirement that the B- and CP-violating interactions must be out of thermal equilib- 
rium In the Standard Model (SM), the sum of B and the lepton number L, B + L, 
is violated anomalously through topologically extended solutions, known as sphalerons. 
In contrast to B + L non-conservation, sphalerons preserve the quantum number B — L. 
The authors in ^ have found that B + L anomalous violation may be large at high temper- 
atures T above the critical temperature of the electroweak phase transition. For values 
of T not much larger than the I^-boson mass, Mw, i.e. T > 200 GeV, up to temperatures 
of T ~ 10^^ GeV, the anomalous B + L rate may exceed the expansion rate of the Uni- 
verse Therefore, any primordial BAU generated at the GUT scale should not rely on 
B + L-violating operators, since sphalerons being in thermal equilibrium will then wash it 
out. The latter appears to be a generic feature of most GUT's, where i? — L- violating terms 
may be suppressed against B + L-violating interactions, thus leading to the net effect of a 
vanishing BAU. On the other hand, it was suggested that the same anomalous B + L 
electroweak interactions may also be utilized to produce the observed excess in B during a 
first-order electroweak phase transition. Given the fact that the experimental lower mass 
bound of the Higgs boson H is about Mh > 80 GeV, this scenario of electroweak baryoge- 
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nesis must now be considered to be rather improbable to explain the observed BAU [P,p!0[ 
within the minimal SM. 

It is therefore important to note that baryogenesis not only provides the strongest 
indication against the completeness of the SM but also poses limits on its possible new- 
physics extensions. An attractive scenario that may lead to a consistent solution to the 
problem of the BAU is the one proposed by Fukugita and Yanagida []TT|,|^ in which the 
baryon number is generated by out-of-equilibrium L-violating decays of heavy Majorana 
neutrinos Ni with masses m^i'^ Tc. Moreover, it was argued |]ll| that the excess in L will 
then be converted into the desired excess in B by means of i? + L- violating sphaleron inter- 
actions, which are in thermal equilibrium above the critical temperature Tc. Many studies 



have been devoted to this mechanism of baryogenesis through leptogenesis |]12|,|13],|14|,|1^,|1| 
over the last years. 

Evidently, possible mechanisms for enhancing CP violation play a decisive role in 



understanding the BAU. In [|IT|, the necessary CP violation in heavy Majorana neutrino 
decays results from the interference between the tree-level graph and the absorptive part 
of the one-loop vertex. Since CP violation originates entirely from the decay amplitude in 
this case, we shall attach the characterization of the e' type to this kind of CP violation, 
thereby making contact with the terminology known from the K^K^ system \Tl\. The 



e'-type CP violation was discussed extensively in the hterature |jT2|0,|Tl]. Provided all 



Yukawa couplings of the Higgs fields to and the ordinary lepton isodoublets are of 
comparable order , baryogenesis through the e'-type mechanism requires very heavy 

Majorana neutrinos with masses not much smaller than 10'' — 10^ GeV. If an hierarchical 
pattern for the Yukawa couplings and the heavy Majorana neutrino masses is assumed 
|l^ , P^ , the above high mass bound may be lifted and the lightest heavy neutrino can have 
a mass as low as 1 TeV. Taking out-of-equilibrium constraints on scatterings involving 
heavy Majorana neutrinos into account |12|, one finds that e' may reach values up to 
10~^ — 10~^ in such a scenario. This must be compared with the usual scenario in [|I^, for 
which e' < 10~^^ for mTVi ~ 1 TeV, and hence very heavy neutrinos are needed to account 
for the BAU. 



In Refs. [p!B|,|5[], the authors pointed out that CP violation in heavy particle decays 
responsible for the baryon (lepton) asymmetry may be further enhanced if one considers 
the absorptive part of the Higgs self-energies, which was though neglected in subsequent 
studies. Since this kind of CP violation may resemble the known mechanism of CP violation 
through K^K^ mixing in the kaon complex |T^, we shall call it hereafter as CP violation 



* The authors in |Q presented an analogous scenario, which was, however, based on anomalous elec- 
troweak decays of exotic Dirac leptons into quarks. 
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of the e type. Exploiting this idea, Botella and Roldan |]I9| gave some estimates for e- 
type CP violation in Higgs decays. They found that the ratio e/e' may indeed be large 
in an extended SU(5) unified model, since the self-energy and vertex contributions to CP 
violation may have different Yukawa coupling structures. Applying this mechanism, the 
authors in |^ concluded that e/e' is of order unity in scenarios for leptogenesis, and 
hence the known results [|ll|] obtained for e'-type CP violation may not need be modified 
drastically. Moreover, using an effective Hamiltonian approach based on the Weisskopf- 
Wigner approximation (WW) , the same authors |^ reached the conclusion that £-type 
CP violation may be rather suppressed, when the two mixed heavy Majorana neutrinos 
are nearly degenerate, and CP violation vanishes completely in the limit, in which the two 
mass eigenvalues of the effective Hamiltonian are exactly equal. 

Recently, there has been renewed interest in the £-type CP violation due to the 
mixing of heavy Majorana neutrinos and the implications of this mechanism for the BAU 
16| . It has been observed in that CP violation can be considerably enhanced through 



the mixing of two nearly degenerate heavy Majorana neutrinos. Using exact solutions for 
the wave functions, which were obtained from diagonalizing the effective Hamiltonian, the 



authors |jT6[ have calculated e and found that it can be larger than e' by two or even 
three orders of magnitude. This result makes the leptogenesis scenario very attractive. 
The enhancement of the CP-violating phenomenon is in agreement with earlier articles on 
resonant CP violation in scatterings involving top quarks, supersymmetric quarks or Higgs 
particles in the intermediate state p2| , p3| , p4| , |25| , as well as with a remark p concerning CP 



violation in the decays of exotic neutral leptons to quarks. 

The existing difference between earlier articles, which found values of e/e' of order 
one and recent authors p^ , |T6i^ , who discovered that e could be even of order 

unity p^ , |25[] , may be attributed to the problem of the proper treatment of two nearly 
degenerate states. It is known that conventional perturbation field theory breaks down in 
the limit of degenerate particles. For example, the wave-function amplitude that describes 
the CP-asymmetric mixing of two heavy Majorana neutrinos, A^i and N2-, say, is inverse 
proportional to the mass splitting m^r^ —nij^^, and it becomes singular if the degeneracy is 
exact [^. Solutions to this problem have been based on the wave-function formalism in 
the WW approximation pT| , p!6[| . Obviously, a more rigorous field-theoretic approach to the 
resonant phenomenon of CP violation through nearly degenerate heavy Majorana neutrinos 
is still necessary. Therefore, it is rather important to provide a field-theoretic solution to 
the problem of e-tjpe CP violation and compare the so-derived results with those found 
with other methods. 



Since the dynamics of e-type CP violation is quite closely related with CP violation 
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induced by particle widths one is therefore compelled to rely on resummation ap- 
proaches, which treat unstable particles in a consistent way. In the context of gauge field 
theories, a gauge-independent resummation approach to resonant transition amplitudes has 
been formulated, which is implemented by the pinch technique p6[. Subsequently, this for- 
malism has been extended to the case of mixing between two intermediate resonant states 
in scattering processes |p4| , p5| . Here, we develop a related formalism for decays, which can 
effectively take into account phenomena of mixing of states during the decay of particles. 

Consequently, our main interest in this paper will be to study the e- and e'-type 
mechanisms of CP violation in some detail, within the framework of an effective field- 
theoretic formalism devised for decay amplitudes. This formalism consistently describes the 
phenomenon of resonantly enhanced CP violation through the mixing of nearly degenerate 
heavy neutrinos and can therefore be applied to any analogous system responsible for 
baryogenesis, in which CP violation is of the e type. The analytic results obtained for e- 
type CP violation with our field-theoretic approach do not display singularities and exhibit 
a physically correct analytic behaviour in transition amplitudes. The fact that resonant 
CP violation through mixing may be of order one l^^^ can lead to scenarios, in which the 
heavy Majorana neutrinos are relatively light with masses as low as 1 TeV. This is roughly 
the highest mass scale, at which the electroweak phase transitions can still occur. Most 
interestingly, this mechanism can produce significant e-type CP violation, even if all Yukawa 
couplings are of the same magnitude and the Majorana masses are of few TeV. This novel 
phenomenological consequence of resonant CP violation has not yet been studied in this 
leptogenesis scenario. Furthermore, it is worth investigating the influence of other possible 
phenomena on this resonant CP-violating mechanism, such as low-energy constraints due 
to the electric dipole moment (EDM) of electron or finite-temperature effects. 

The paper is organized as follows. In Section 2, we describe minimally extended 
models that include heavy Majorana neutrinos. At low energies, these models amount to 
adding isosinglet neutrino states to the field content of an effective one-Higgs doublet model. 
Such scenarios may be embedded into certain SO(IO) and/or Eg unified theories, which 
can naturally predict nearly degenerate heavy Majorana neutrinos as light as 100 GeV. 
Moreover, we discuss the renormalizability of the effective model. In Section 3, we present 
a resummation formalism for resonant transitions between fermions in decay amplitudes. 
Furthermore, we illustrate some of the advantages of this approach, when compared to 
other existing methods. Making use of our formalism, we calculate the analytic expressions 
of the relevant transition amplitudes in Section 4. In Section 5, we give estimates of possible 
constraints coming from low-energy data, such as the EDM of electron, which turn out to 
be quite weak in order to rule out our leptogenesis scenario. In Section 6, we present the 
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Boltzmann equations relevant for the evolution of the leptonic asymmetry in the effective 
model, and give numerical estimates and comparisons for the BAU generated via e- and/or 
e'-type CP violation. We also discuss the implications of finite temperature effects for the 
resonant phenomenon of CP violation and find that such a phenomenon can still be viable. 
We draw our conclusions in Section 7. 



2 Heavy Majorana neutrino models 

Heavy Majorana neutrinos may naturally be realized in certain GUT's, such as SO (10) 
^ and/or Eg |]2^ models. Nevertheless, these models will also predict several other 



particles, e.g., leptoquarks, additional charged and neutral gauge bosons {W^ and Zr), 
which may deplete the number density of heavy neutrinos Ni through processes of the type 
NiCpi W^* uudu and so render the whole analysis very involved. If these particles 
are sufficiently heavier than the lightest heavy Majorana neutrino and/or the temperature 
of the universe this problem may be completely avoided. Since we wish to simplify 



our analysis without sacrificing any of the essential features involved in the study of the 
BAU, we shall consider a minimal model with isosinglet neutrinos, which is invariant under 
the SM gauge group, SU(2)/^(8)U(l)y. Then, we will present the relevant Lagrangians that 
govern the interactions of the heavy Majorana neutrinos with the Higgs fields and the 
ordinary leptons. Also, we will identify the non-trivial CP-violating phases of the model 
and pay special attention to the one-loop renormalization of the Yukawa couplings. 



As has been mentioned above, certain SO(IO) p7| , |28| and/or Eg [29| models natu- 
rally predict the existence of heavy Majorana neutrinos. In SO(IO) models, an attractive 
breaking pattern down to the SM may be given schematically in the following way: 

SO(IO) ^ G422 = SU(4)ps ® SU(2)fl ® SU(2)z. 

G3221 = SU(3)e ® SU(2)r ® SU(2)z. ® U{l)iB-L) 
^ SM = G32i = SU(3)®SU(2)L®U(l)y, (2.1) 



where the subscript PS refers to the Pati-Salam gauge group [^. The spinor representation 
of SO(IO) is 16 dimensional and its decomposition under G422 is given by 

ai22 : 16 ^ (4,1,2) © (4,2,1). (2.2) 



As can be seen from Eq. ( |2.2| ), it is evident that SO(IO) can accommodate right-handed neu- 
trinos, since it contains the left-right symmetric gauge group SU(2)/j(S)SU(2)l®U(1)(b_l). 
There are several Higgs-boson representation that can give rise to the breakdown of G422 
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and G3221 down to the SM gauge group 0^21 p8|j3T|] . In Ee models the 27 spinor 
representation decomposes into 16 © 10 © 1 under SO(IO), which leads to four singlet 
neutrinos per SM family: one neutrino as isodoublet member in 16, two neutrinos as 
isodoublet members in 10, and one singlet neutrino in 1. In these models, two of the four 
isosinglets can have Majorana masses of few TeV [^, depending on the representation of 



the Eg Higgs multiplets, whereas the other two are very heavy with masses of the order 
of the unification scale. Possible flavour structures for the isosinglet neutrino mass matrix 
resulting from the above two representative unified models will be discussed below. 

The minimal model under consideration extends the SM field content of the three 
lepton and quark families by adding a number ur right-handed neutrinos with i = 
1,2,..., Tin. Even though in Eg models the active isosinglet neutrinos may be more than 
three, in the SO(IO) models mentioned above the symmetric case of having one right- 
handed neutrino per family turns out to be quite natural, i.e., rtn = 3. Therefore, we shall 
not specify the number of z/^^j in the following. To be specific, the leptonic sector of our 
minimal model consists of the fields: 

with / = e,ix,T. Since for temperatures T ^ v, one has v{T) = 0, where v(T) is 

the vacuum expectation value (VEV) of the SM Higgs doublet $ at temperature T (with 
V = f (0)), the only admissible mass terms are those of the Majorana type and are given 
by the Lagrangian 

-Ui = ^ E (^^.M^.z/H, + (2.3) 

^ i,j=l 

Here, the superscript C denotes the operation of charge conjugation, which acts on the 
four- component chiral spinors ipL and ipR as (V'l)'" = PrCiP^ and (V'i?)'" = PlC^'^ , where 
Pl{r) = [1 ~ (+)75]/2 is the chirality projection operator. In Eq. (^]3|), is a n^j x tir 




dimensional symmetric matrix, which is in general complex. The isosinglet mass matrix M'^ 
can be diagonalized by the unitary transformation U'^M'^U = , where U is a ur x 
dimensional unitary matrix and M'^ is a positive diagonal matrix containing the ur heavy 
Majorana masses. Furthermore, the ur mass eigenstates Ni are related to the flavour 
states URi through URi = PRYJj=iUijNj. In the basis, in which the isosinglet neutrino 
mass matrix is diagonal and equals , the Yukawa sector describing the interaction of 
the heavy neutrinos with the Higgs doublet and the ordinary leptons reads 

>Cy = - EE h) ( ~ "^'^'"^ \ N, + H.c. (2.4) 

1=1 j=i \ -X 
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The CP-even Higgs field if, the CP-odd Higgs scalar x° and the charged Higgs scalars 
given in (|2.4|) are all massless at high temperatures. In the limit T ^ 0, H represents the 
massive SM Higgs boson, whereas and ^ire the massless would-be Goldstone bosons 
eaten by the longitudinal degrees of freedom of the gauge bosons Z and W^, respectively. 
If the SM Higgs field H is very heavy with mass of order 1 TeV at T = and also close to 
the mass of the decaying heavy neutrinos, then mass Higgs effects may not be negligible. 
Therefore, we shall initially keep the full Mh dependence in our calculations and then 
present analytic results for the limiting case Mh = 0. 

Before counting all the non-trivial CP-violating phases for the most general case, it 
may be more instructive to discuss first some simple models that could predict degenerate 
or almost degenerate heavy Majorana neutrinos and CP violation. To this end, we consider 
a minimal scenario, in which the fermionic matter of the SM is extended by adding two 
right-handed neutrinos per family, e.g., vir and {Sil)^ with I = e,fi,T. Such a scenario 
may be derived from certain SO(IO) and/or Eg |^ models. For our illustrations, we 



neglect possible inter-family mixings. Imposing lepton-number conservation on the model 
gives rise to the Lagrangian 




-C = -{Sl, (i^Rr) . \ + hRii^L, D^i^B^ + H.C., (2.5) 



where $ = ia2^ is the isospin conjugate Higgs doublet and (J2 is the usual Pauli matrix. 
Even though the mass and coupling parameters may be complex in such a scenario, the 
phase redefinitions of the fields, 

lyL-^e^'^-^UL, h^e^^HL, ur-^c'^^ur, Sl^c^^^Sl, (2.6) 

can, however, make them all real. This model is CP invariant, unless one allows for a non- 
trivial mixing among generations Moreover, the model preserves the lepton number 



and hence cannot produce any excess in L through heavy neutrino decays. 

There are two equivalent ways to break the L and CP invariance of the Lagrangian 
in Eq. ( |2.5| ): One has to either (i) introduce two complex L- violating mass terms of the 
kind f^Ri^Ri^R and fiiS^SL, where both fiR and /i^ are complex, or equivalently (ii) add the 
L-violating coupling hR{vL, Il)^{Sl)^ and include the L-violating mass parameter, e.g., 
f^Ri^R'^R- These are the minimal enlargements that can assure L and CP violation on the 
same footing in this simple two-isosinglet neutrino model. In fact, the necessary conditions 
for CP invariance in these two scenarios are written down 

(i) \hR\'Qm{M*^iJLm) = 0, 

(ii) '^m{hLh*J,^IRM*) = 0. (2.7) 
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It can easily be checked that the two equahties in Eq. ( p.7| ) are invariant under the phase 
redefinitions of the fields given in Eq. ( p.6|) . We must remark that the L-breaking pa- 
rameters fiL and fiR are generally much smaller than M within Eg scenarios. The ori- 
gin of these parameters are usually due to residual effects of high- dimensional operators 
involving super-heavy neutrinos [^. The typical size of the L- violating parameters is 
fiL, fiR ~ M'^/Mpianck, M'^/Mx OT M'^/Ms, whcrc Ms ^ 10^^ is some interme- 
diate see-saw scale. As a consequence, such effective minimal models derived from Eg 
theories can naturally predict small mass splittings for the heavy neutrinos A^^i and N2. 
To a good approximation, this small mass difference may be determined from the pa- 
rameter xjy = 'mjq.^/mN^ _ 1 ^ /Ul/M or jiji/M. For instance, if M = 10 TeV and 
/^L = /^i?, = /Mxi one then finds xn ~ lO^"*^^ — 10~^^. As we will see in Section 4, these 
small values for the mass difference xn can produce large CP asymmetries in the heavy 
neutrino decays. 

In order to deduce the sufficient and necessary conditions for the most general struc- 
ture of the two right-handed neutrino model, one must consider the systematic approach 
presented first in [Q, which is slightly different from the procedure outlined above. In 



this approach [ p^j33[l , one looks for all possible weak-basis independent combinations that 



can be formed by Yukawa couplings and the neutrino mass matrix Af, and are simultane- 
ously invariant under generalized CP transformations of the fields. These generalized CP 
transformations of the fermion fields may include unitary fiavour rotations, apart from the 
phase redefinitions mentioned above. Further details may be found in Ref. Thus, for 
the model at hand, we find that the sufficient and necessary condition for CP invariance is 

QmTT{h^hM''^M''M''^h^h*M'') = mN^mN^{m%^ - mlj '^mihah^^Y = 0. (2.8) 



The 1x2 dimensional matrix h in Eq. ( p. 81 ) contains the Higgs Yukawa couplings, which 
are defined as hij in Eq. ( p.4|) , i.e., in the physical mass basis where = . One 
can show that Eq. (|2.8|) is consistent with the conditions in Eq. ( p.7|) for the two special 
cases discussed above. From Eq. (p^), one readily sees that only one physical CP- violating 
combination is possible in this minimal model and CP invariance is restored if m^r^ = mjsf^ 
provided none of the isosinglet neutrinos is massless. The above considerations may be 
extended to models with more than two right-handed neutrinos and more than one lepton 
families. In this case, there may be more conditions analogous to Eq. ( p.8|) , which involve 
high order terms in the Yukawa-coupling matrix h. However, not all of the conditions are 
sufficient and necessary for CP invariance. Instead of undertaking the rather difficult task 
to derive all CP-invariant conditions, we note in passing that the total number N'cp of 
all non-trivial CP- violating phases in a model with ni weak isodoublets and ur neutral 
isosinglets is Afcp = tiLijiR — 1) 
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Fig. 1: One- loop graphs contributing to the renormahzation of the couphngs x l^i^ 
X^viNi and HuiNi. 

Since CP violation in Ni decays will necessitate non-zero one-loop absorptive parts of 
vertex and self-energy graphs as will be seen in Section 4, one should also have to address 
the issue of renormahzation of the dispersive counterparts (see also Fig. 1). In order to 
check explicitly that our minimal model leads indeed to consistent renormalizable results, 
we shall adopt the following strategy in our analysis. First, we fix the renormahzation of 
all Higgs Yukawa couplings hij from the decay mode Ni l^X~ ■ this way, we determine 
the counter-terms (CT) of /i/j, 5hij. Then, we show that all ultra-violet (UV) divergences 
cancel in the partial decays uix^ and A^'j uiH. For this purpose, we first express 

all bare quantities in terms of renormalized ones as follows: 



i'=i ^ 



'■L 



l'=l ^ 

(l + ^SZ^)^, hi = hi, + Shi, 



(2.9) 



The superscript '0' in Eq. ( p.9|) indicates that the field or coupling parameter is unrenor- 
malized, whereas quantities without this superscript are considered to be renormalized. In 
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addition, the CT 6Z^ collectively denotes the wave-function renormalization constants of 
all components of the Higgs doublet $ (or $), i.e., the fields x^, X° ^ind H. The divergent 
part of all the Higgs wave-function renormalizations, 6Z^^^, has been found to be universal. 
Expressions showing the universality of together with other relevant one-loop analytic 
results are relegated to Appendix A. Taking the relations in Eq. (|2.9|) into account, we find 
that the renormalized Lagrangian in ( p. 41) gets shifted by an amount 

T riL UR r, riL . ^ 

-^^y = o EE + '^^^ + T^^Zl^i + j:^Z^k)U^N, + H.c, (2.10) 

^ 1=1 j=l '^Ij V=l k=l 

where Li = (z//^, 1^)'^ and 6Z'^ = {5Z\5Z^). From Fig. 1(a), it is easy to see that the 
one-loop correction to the coupling x^^^ can only occur via a AL = 2 Majorana mass 
insertion, owing to charge conservation on the vertices. Naive power counting may then 
convince oneself that the one-loop irreducible vertex x^A^/ is UV finite. Less obvious is the 
UV finiteness for the proper couplings x°A^z/ and HNu, which is shown in Appendix A. 

As has been discussed above, we shall now determine the Yukawa coupling CT's 6hij 
from the renormalization of the coupling x^^^- Requiring that all UV terms are absorbed 
into the definition of hij, we obtain 

1 "L "fl 

= - - ( hi,5Z^- + hi'jSZl* + Y: hik^ZZ ) . (2.11) 
^ v=i k=l 

We observe that 5hij may be separated into two terms: The wave-function term 5Z-^-, 
which is flavour independent, and the rest, which depends on the flavour and the wave- 
function CT's 5Z\i, and ^Zf^ . If we had renormalized the Higgs Yukawa couplings from 
the decays uH, the only difference in Eq. (|2.11|) would have been the appearance 

of the CT's 6Zh and SZ^, in place of SZ^~ and SZ^,, respectively. Also, for the decay 
A^ z/x", one has to include the wave-function renormalization of ^^x°^ instead of 
6Zh in the decay A^ —>■ vH . Therefore, consistency of Yukawa coupling renormalization 
requires that differences of the kind 5Z\i, — 5^;*^/, 5ZyO — 6Zh and SZ^- — 5Zh must be 
UV safe. In Appendix A, it is shown that all these CT differences are indeed UV finite 
and vanish in the limit of Mh 0. This completes our discussion concerning the one- 
loop renormalization of heavy Majorana neutrino decays. In the next section, we shall 
explicitly demonstrate how the renormalization presented here gets implemented within 
our resummation formalism for unstable particle mixing. 
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3 Resummation approach for two-fermion mixing 



If a Lagrangian contains unstable particles, then these fields cannot be described by free 
plane waves at times t ±00 and hence cannot formally appear as asymptotic states in 
the conventional perturbation field theory. Within a simple scalar theory with one unstable 
particle, Veltman showed that, even if one removes the unstable particle from the initial 
and final states and substitutes it in terms of asymptotic states, the so-truncated S-matrix 
theory will still maintain the field-theoretic properties of unitarity and causality. Our main 
concern in this section will be to present an approach to decay amplitudes that describes 
the dynamics of unstable particle mixing. Hence, in such a formulation, finite width effects 
in the mixing and decay of non-asymptotic states must be taken into account. This will 
be done in an effective manner, such that the decay amplitude derived with this method 
can be embedded in an equivalent form to a transition element p6| , p5| in agreement with 
Veltman's S-matrix approach. This effective field-theoretic approach is equivalent to that 
of the decay of an initial pure state, such as the states or K'^, which is initially produced 
by some asymptotic states in kaon experiments, e.g., in piT~ or pp collisions [0. Since 



the time evolution of the decaying system is effectively integrated out over all times, the 
resummed decay amplitudes derived with this field-theoretic method will not display any 
explicit time dependence. 

The discussion in this section is organized as follows. First, we briefiy review the 
theoretical description of the mixing between stable particles in a simple scalar theory 
within the framework of the Lehmann-Symanzik-Zimmermann formalism (LSZ) After 
gaining some insight, we extend our considerations to the mixing between two unstable 
scalars. The effective field-theoretic method developed for the scalar case can then carry 
over to the case of mixing of two unstable fermions, with the help of which e-type CP 
violation will be calculated in Section 4. 

Let us now consider a field theory with N real scalars 5^°, with i = 1,2, . . . , N. We 
shall assume that the scalars are stable to a good approximation and neglect possible finite 
width effects of the particles. The bare (unrenormalized) fields and their respective 
masses may then be expressed in terms of renormalized fields Si and masses Mj in the 
following way: 

= Zfs, = (5,, + ISZ,,)S,, (3.1) 
(M°)2 = Mf + 6M^ . (3.2) 

Here and in the following, summation is understood over repeated indices that do not 
appear on both sides of an equation. In Eqs. (|3.1|) and (|3.2| ), Z^Y^ and 6Mi are the wave- 
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function and mass renormalization constants, respectively, which can be determined from 
renormahzation conditions imposed on the two-point correlation functions, Uijlp^), for the 
transitions 5*^ Si in some physical scheme, such as the on-mass-shell (OS) renormaliza- 



tion scheme 37 



Si 



lim 




Fig. 2: Diagrammatic representation of the renormalized n — 1-non-amputated amplitude, 
Sj^,.., and the LSZ reduction formalism. 

It will prove useful for the discussion that follows to give the relation of the pole parts 
between the unrenormalized scalar propagators Ajj(p^) and the renormalized ones Ajj(p^). 
The two pole parts are related through p7| 



Using the LSZ formalism shown schematically in Fig. 2, one can deduce the renormalized 
— 1-non-amputated amplitude, Si^,,,, for a fixed given external line i, from the correspond- 
ing unrenormalized n-point Green function Gj,,., where n is the total number of external 
lines. In this way, we have 

5,... = hm G,..ZJi"'\p^ - Ml) 

p2^Mf ' — M^ ■' ■' 

= hm r-^Zli' , (3.4) 

where T^"^^ denotes the amplitude amputated at the k external leg. Clearly, the LSZ re- 
duction procedure outlined above can be generalized to all external legs, thus leading to the 
physical (renormalized) S'-matrix element 5'j^.,,j^, which governs the transition amplitude 
of n asymptotic states. 

Let us now consider the mixing of two neutral unstable scalars p^ , |25| , e.g., Si and 
5*2. Since we are interested in studying the width effects of these particles, we have first 
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to calculate all the SiSj Green functions, with i,j = 1,2. After summing up a geometric 
series of the self-energies Ilij{p^), the full propagators may be obtained by inverting the 
following inverse propagator matrix: 



(3.5) 



The result of inverting the matrix in Eq. ( p.5|) may be given by 







P'- 




)^ + nn(/)- 


A22(/) 




p'- 


(M° 


)' + n22(/)- 


Al2(/) 


= A2l(/ 


) = 


-ni2(s) 



n?2(P^^ 



p2_(M0)2 + n22(p2) 

p2-(M0)2 + nn(p2) 



(3.6) 
(3.7) 



1 ~i 



n?2(P^^ 



(3.^ 



where Ili2{p'^) = 1121 (p^). Moreover, we find the useful factorization property for the off- 



diagonal {i 7^ j) resummed scalar propagators 
A,,(p2) = -A,,(p2) 



p2 - (MO) 2 + %(p2 



p2 _ (^0)2 + n,,(p2) 



A,,(p2) . (3.9) 



The resummed unrenormalized scalar propagators Aij{p'^) are related to the respective 
renormalized ones Ajj(p2) through the expression 



(3.10) 



where Aij^p"^) may be obtained from Eqs. ( p.6|) -(3.8), just by replacing AfO with Mi and 
Hijlp"^) with IIjj(p2). Note that the property given in Eq. (|3.9| ) will also hold true for the 
renormalized scalar propagators Aij{p'^). Taking expressions ( |3.9| ) and ( ^.10[ ) into account, 
we can derive the resummed and renormalized transition amplitude, denoted here as S'j,,,, 
for the external leg i which now represents an unstable particle. This can be accomplished 
in a way analogous to Eq. ( |3.4| ), viz. 

= lnn^^rr'Zl!:AUp')Zlfz;^/'^AT,\p') 



lim 



T' 



amp /7 1 /2 



ki 



^ampr^l/2 nmj(p2)(l — (5^.) 



km 



p-^-Ml + Ummip^ 



Mf-M] + %,{M! 



(3.11) 
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where Si,„ and Sj,,, are the renormahzed transition elements evaluated from Eq. ( p.4| ) in 
the stable-particle approximation. One should bear in mind that the OS renormalized 
self-energies Ilji{Mf) in Eq. ( 3.11| ) have no vanishing absorptive parts, as renormalization 
can only modify the dispersive (real) part of these self-energies. The reason is that the CT 
Lagrangian must be Hermitian as opposed to the absorptive parts which are anti-Hermitian. 
In fact, these additional width mixing effects are those which we wish to include in our 
formalism for decay amplitudes and are absent in the conventional perturbation theory. It 
is also important to observe that our approach to decays is not singular, i.e., Si,„ displays 
an analytic behaviour in the degenerate limit — ^ Mj, because of the appearance of the 
imaginary term iQmIljj{Mf) in the denominator of the mixing factor present in the last 
equality of Eq. ( p.ll ). Finally, we must stress that the inclusion of these phenomena has 
been performed in an effective manner. Since the decaying unstable particle cannot appear 
in the initial state [^, the resummed decay amplitude must be regarded as being a part 
which can effectively be embedded into a resummed ^-matrix element |2^. This resummed 
S'-matrix element describes the dynamics of the very same unstable particle, which is 
produced by some asymptotic states, resides in the intermediate state, and subsequently 
decays either directly or indirectly, through mixing, into the observed final states. 

It is now straightforward to extend our considerations to the case of mixing between 
two unstable fermions. Following a line of arguments similar to those presented above, we 
consider a system with two unstable fermions, call them /i and /2. As usual, we express the 
bare left- and right-handed chiral fields, f^^ and (with i = 1, 2), in terms of renormalized 
fields as follows: 



Jl 



Li 



^Lij JLj 





Ri 



Rij JRj ' 



where Z]/-^ is the wave-function renormalization constant for the left 



(3.12) 
right-) handed 

chiral fields, which may be determined from the fermionic self-energy transitions fj — ^ 
/j, Sjj(^, e.g., in the OS renormalization scheme [^. Analogously with Eq. ( p.5|) , the 
resummed fermion propagator matrix may be obtained from 



S; 



^-771? + Sll(^) 



^211 



221 



(3.13) 



where 2 are the bare fermion masses, which can be decomposed into the OS renormalized 
masses mi^2 and the CT mass terms Smi^2 as = ^1,2 + ^^1,2- Inverting the matrix- 
valued 2x2 matrix in Eq. (|3.13| ) yields 



5*22 (^) 



221 



^21 



^211 



- + S22(^) 



^- + Sii(^) 



S2l(^) 



-1 



(3.14) 
(3.15) 
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S2M 



rrio 



^ - m? + Sn(^)] Si2(:^)522(^) 
5'22(j^) S2i(^) ^ - + Slid 



-1 



- + S22(^) S2i(;^)5ii(:^) . 



(3.16) 



(3.17) 



From Eqs. ( |3.16 ) and ( 3.17 ), it is now easy to see that the resummed propagator matrix is 
endowed with a factorization property analogous to Eq. ( p.9|) . There is also an analogous 
connection between the renormalized and unrenormalized resummed propagators, which 
may be cast into the form 



7l/2t p N 



(3.18) 



where the caret on Sij{'f) refers to the fact that the resummed fermionic propagators have 
been OS renormalized. By analogy, the renormalized propagators Sij{]f) may be recovered 



from Sij{p) in Eqs. ( |3.14| )-( |3.17| ), if one makes the obvious replacements: 



rrii and 



Employing the LSZ reduction formalism, one can derive the resummed decay ampli- 
tude. Si, of the unstable fermion fi —>■ X , in a. way similar to what has been done for the 
scalar case. More explicitly, we have 



SiUi{p) 



rj-amp 
^k... 



{^]!km Pl 
l/2t 



7I/2 p 
^Rkm ^R) "^mn 



AzirpR + z 



SiUi{p) 



-l/2t 
Rji 



Pt 



y\l2\ p yll2\ 



'Rnj 



Pl) 



(1 - bij)Sj f^ji^f) 



-1 



u,{p). (3.19) 



Again, Si represent the respective renormalized transition amplitudes evaluated in the 
stable-particle approximation. The amplitudes Si also include all high n-point functions, 
such as vertex corrections. On the basis of the formalism presented here, we shall calculate 
the CP asymmetries in the decays of heavy Majorana neutrinos in Section 4. 

Finally, we wish to offer a comment on other approaches, which are used to analyze 
the phenomenon of CP violation through particle mixing [^,|l^. Recently, this issue has 
been studied in , within a S-matrix amplitude formalism related to the one discussed in 
this section. Obviously, the approach based on the diagonalization of the effective Hamil- 



tonian piipi)| is very helpful to describe e-type CP violation, if the effective Hamiltonian 
is diagonalizable through a similarity transformation, as is the case for the known K^K^ 
system. However, if the effective Hamiltonian has mathematically the Jordan form, when 
expressed in a K^K^-\\ke basis, then it can be shown to be non-diagonalizable via a simi- 
larity transformation. In this case, the complex mass eigenvalues of the two mixed non-free 
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particles are exactly equal and, most importantly, CP violation through particle mixing 
reaches its maximum attainable value |25[ . 

To give a specific example, let us consider the following effective Hamiltonian for the 
mixing- system of two nearly degenerate heavy neutrinos A^^i and N2: 



mi - 



7712 - ^22(1 



m]\r + a — i\b\ 
-ib* 



-ib 



i\b\ 



(3.20) 



in the approximation ^ ^ rriN ^ nii ^ 1112. In Eq. (|3.20|) , the parameters a and b are 
real and complex, respectively, and mi = niN + a, 1112 = rriN — a. The complex parameter 
b represents the absorptive part of the one- loop neutrino transitions Ni Nj. Unitarity 
requires that the determinant of the absorptive part of 7Y(je() be non- negative. For the 
effective Hamiltonian (|3.20|) , the corresponding determinant is zero. Such an absorptive ef- 
fective Hamiltonian naturally arises in the one-lepton doublet model with two right-handed 
neutrinos. In the limit a — >■ the two complex mass eigenvalues of 7i(^) are exactly de- 
generate and equal to m^r — As has been shown in [^, the effective Hamiltonian 
cannot be diagonalized by a non-unitary similarity transformation in this limit, i.e., the 
respective diagonalization matrices become singular. 

Since our effective field-theoretic method does not involve diagonalization of the ef- 
fective Hamiltonian through a similarity transformation, such singular situations are com- 
pletely avoided, leading to well-defined analytic expressions for the decay amplitudes. Fur- 
thermore, whenever referring to heavy Majorana neutrino masses in the following, we shall 
always imply the OS renormalized masses within the conventional perturbation field theory, 
which differ from the real parts of the complex mass eigenvalues of the effective Hamilto- 
nian. In the presence of a large particle mixing, the corresponding eigenstates of the latter 
are generally non-unitary among themselves, whereas the eigenvectors of the former form 
a well-defined unitary basis, upon which perturbation theory is based. In this context, it 
is important to note that these field-theoretic OS renormalized masses are those that enter 
the condition of CP invariance given in Eq. (|2.8|) . Therefore, we can conclude that, if the 
two complex mass eigenvalues of the effective Hamiltonian for the mixed heavy Majorana 
neutrinos are equal, this does not necessarily entail an equality between their respective 
OS renormalized masses, and, hence, absence of CP violation as well p5|. 



4 CP asymmetries 

It is now interesting to see how the resummation formalism presented in Section 3 is applied 
to describe e-tjpe CP violation in heavy Majorana neutrino decays. The same formalism 
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can be used for the inverse decays, which occur in the formulation of the Boltzmann equa- 
tions (see also Section 6). For completeness, we shall include e'-type CP violation in our 
analysis, which originates entirely from the one-loop ^IN irreducible vertex, and display 
plots with numerical comparisons between the two kinds of CP-violating contributions 
mentioned above. Moreover, we wish to address briefly the issue pertaining to the problem 
of flavour-basis invariance of the CP asymmetries. 

^.^ <|)t 




L 



Fig. 3: e- and e'-type CP violation in the decays of heavy Majorana neutrinos. 



Let us consider the decay A^^i — > l~x^ ^ model with two-right handed neutrinos, 
shown in Fig. 3. The inclusion of all other decay modes will then be straightforward. To 
make our resummation formalism more explicit, we shall first write down the transition 
amplitude responsible for e-type CP violation, denoted as Tlf \ and then take CP- violating 
vertex corrections into account. Applying Eq. ( p.l9| ) to the heavy neutrino decays, we have 



T 

J- 1\ 



Ni 



hiiUiPnUN, - ihnUiPn ^ - mTv^ + iS^/(^) 



In Eq. (4J^), the absorptive part of the one- loop transitions Nj ■ 
the general form 

= A,,(/) + a:^{p') ^Pn 

where 

'^^^ ' 327r L2 2 V 
In the limit Mh —>■ 0, one finds the known result |lT5|,p!6| A. 



^'lltiMN, . (4.1) 
Ni, with i,j = 1,2, has 



MfJ^2 



(4.2) 



(4.3) 



hiiihlJilQ-K). On the other 



hand, the CP-transform resummed amplitude describing the decay A'^i 
written down 



Nil 



IS 



^abs , 



= uiPlUni - ih*i2 uiPl [ ^ - tun^ + i^lTi:^) 
where the charge-conjugate absorptive self-energy is given by 



PlVi 



^21 ( 



)UNi 



Mp') ^Pr + AUp") ^Pl 



(4.4) 



(4.5) 
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In deriving the last step of Eq. ( [4.4| ), we have made use of the known identities: u{p, s) = 
Cv^{p,s) and C'-y^C^^ = — 7j- The expressions in Eqs. (^4.1|) and (|4.4| ) may be simphfied 
even further, if the Dirac equation of motion is employed for the external spinors. Then, 
the two resummed decay amplitudes, TJ^^'* and T^^i^ take the simple form 



T^l^ = uiPrUnA hii - ihi2 „ ■ a'' \'? ""^""^ , (4.6) 



(1 + iA22)A*^ + mN^mN 2 Mi 
m\{l + iA22Y - rn%^ 
m% (1 + iA22)A2i + mN^niN^Al 



T„, = u,P,u,,[h, - ,h,, ■ (4.7) 



The two CP-conjugate matrix elements differ from one another in having complex conjugate 
Yukawa couplings to each other and scalar currents with opposite chirality. Furthermore, 



the respective transition amplitudes involving the decays N2 I x'^i '^nJ j ^ind N2 — > I'^x' 



T^nI^ be obtained from Eqs. ( |4.6| ) and (|4.7|), just by interchanging the indices '1' and 
'2' everywhere in the above two formulas. 

We shall now focus our attention on studying the e- and e'-type mechanisms of CP 
violation in heavy Majorana neutrino decays. For this purpose, we define the following 
CP-violating parameters: 

i'r(^)|2 |T=(^)|2 

\-^Ni\\±Nj 



Sn. = ^ ^g^, for ^ = 1,2, (4.^ 



= ''^^ ' + ''^^ ' (4.9) 

Correspondingly, one can define the CP- violating parameters e'j^^ and e^, which may quan- 
tify CP violation coming exclusively from the one-loop irreducible vertices. In Eqs. ([4.8|) 
and ([4. 91), the parameters ^at. and En share the common property that do not depend on 
the final state that Ni decays, despite the fact that the individual squared matrix elements 
do. In general, both e- and e'-type contributions may not be directly distinguishable in the 
decay widths r(A^i — > /^x^)? unless one has s'^. and vice versa, for some range of 

the kinematic parameters. Evidently, the physical CP-violating observables we are mainly 
interested in are 

^ r(iv. ^ L#t) - r(iv. ^ Ly 



- ^^2 — WVt , r^i\ , v-2 



(4.11) 



Etir(iv. ^ i.$t) + Etir(iv. ^ i.^$) 

where L refers to all fermionic degrees of freedom of the leptonic iso doublet that heavy 
Majorana neutrinos can decay. Nevertheless, the parameters ^at., En and e'j^^ defined above 
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are helpful to better appreciate the significance of the two different mechanisms of CP 
violation. 



To elucidate our resummation formalism further, it may be useful to calculate the 
analytic form of the parameters ^at. for the interesting case of nearly degenerate heavy 



Majorana neutrinos. Therefore, we shall consider the approximations Amj^ = mj 



m 



^ "^ATi ~ ""^ATj define the parameter r^r = Am% / (171^-^171^^) . Moreover, we neglect 
high order Yukawa couplings of 0{hfj) at the amplitude level. It is then not difficult to 
find the approximate expressions 

_ ^m{h*i^hi2hl-yhv2) , . 

^m{h*i^hi2h*i,^hv2) , . 

" — Mi^^ — WTW.- 

The difference between the above expressions and those obtained within the framework of 
the ordinary perturbation theory is that the regulating terms Af^ and A22 are absent in 
the latter. Clearly, if these finite width terms were not consistently taken into account, 
this would cause a singular behaviour when the degeneracy between the two heavy Majo- 
rana neutrinos is exact |T^,T^,TB|. On physical grounds, however, this should not be very 



surprising, since the only natural parameter that can regulate such a singularity is the 
finite width of the heavy Majorana neutrinos. Therefore, one of the main advantages of 
our approach is that the dynamics of CP violation through heavy-neutrino mixing can be 
properly described by giving rise to physically well-behaved analytic expressions. 

Another important point, also reported in is the fact that both parameters e^i 



and contribute constructively to CP violation. Technically speaking, this can be seen 
as follows. As has been mentioned above, the expression in Eq. ( 4.1^ ) may be obtained 



from Eq. ([4.12|) , if one replaces the heavy-neutrino index '1' with '2' everywhere in that 



formula. As a result, the CP- violating combination of Yukawa couplings, Qm{hi^hi2hJ,ihi'2) , 
fiips sign, which gets compensated by a similar sign fiip in the parameter r^. Another 
significant feature of our analytic results in Eqs. (|4.12|) and (|4.13|) is that both and 



vanish in the limit rriN-^ — > m^rj, which is consistent with the requirement of CP invariance 



given in Eq. ( |2.8| ). We have checked that this property of CP invariance persists even if 
one calculates Eni and En2 exactly from the expressions (|4.6| ) and (|4.7| ). In addition, we 
have verified that En2 {^Ni) vanishes in the limit rriNi (rriN^) 0, as is also prescribed by 
Eq. ( |2.8| ). As a consequence, our analytic expressions for £-type CP violation are indeed 
proportional to the fiavour-basis invariant combination of CP non- invariance (c/. Eq. (|2.8| )), 
mNimN2{m%^ - m%^)'^m{h*i^hi2h}^hv2). 
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In order to make our analysis complete, it is important to include the contributions 
from e'-type CP violation, since they may be significant for very large differences of heavy 
neutrino masses, e.g., for niN^ — niN^ ~ tjin-^ or m^^. In this regime, the e-type terms 
are suppressed by the large virtuality of the heavy neutrino propagator and so become 
comparable to the e'-type terms It is now useful to define the function 



Fix, a) = \fx 



-«-(l + x)ln(^ — ^^-^^^ 

If one sets a = in Eq. ( [4.14| ), then F[x,a) reduces to the known function j[x) = 
\/x[l — [1 + x) ln(l + found in [0. Here, we are only interested in the L- violating 

absorptive parts of the one-loop vertices x^^^ii X^^i^i ^^id HuiNi, shown in Figs. l(a)-(c). 
The complete analytic expressions are calculated in Appendix A. For later convenience, we 
also assume that the external decaying heavy Majorana neutrinos are off-shell. Having this 
in mind, we find the off-shell absorptive couplings 

V?»,W = -^1^ ^P.f(^.o). (4.15) 



~ 327rv^ 



^ftK^.o) (4.16) 



To make contact with the e'j^_ expressions existing in the literature ||TT],|T2|,|l3|,|T4l , we compute 
the e'-type CP asymmetry in the conventional perturbation theory, using Eqs. ( |4.15| ) and 
( [4.16| ) and neglecting wave-function contributions. Considering all decay channels for the 
decaying heavy Majorana neutrino, e.g., Ni, we obtain 



+ j(l-%W^,0) (4.17) 



In the vanishing limit of the Higgs-boson mass, the above formula simplifies to [|rT],|T^|T^,|I^ 



= fiz::^) ■ (4-18) 



Vl 



Unlike e^^^, e'^^ does not vanish in the degenerate limit of the two heavy Majorana neutrinos 
Ni and N2. However, when the value of m^Vi approaches to that of mjq^, the e'-type part 
of the transition amplitude squared for the A^^i decay becomes equal but opposite in sign 
with the respective one of the N2 decay. Thus, these two e'-type terms cancel one another, 
leading to a vanishing result for the CP- violating parameter e'j^ (7^ e'^^ + s'j^^), which may 
be defined analogously to Eq. 
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We are now in a position to implement the vertex corrections of the Yukawa couphngs 

('^) ^r^A T^^^ rr;,r^v^ ^ "p^c (\A M\ c^^A (\A A\\ fesp ect ivcly. 



to the resummed amphtudes, T^^ and T^^, given in Eqs. ( ^TTP and 
Taking Eqs. ( |4.15| ) and ( 4.1(jD into account, we find 



T 



hi2 + iV) 



abs I 
12 \ 



(4.19) 
(4.20) 



Ul 



where we have simphfied the notation of the off-shell one-loop vertices as Vjf'^i^)- The 
vertex functions V"j^*(^) are the charge conjugates of Vjf'^^'^) and may hence be recovered 
from Eqs. ( [4.15[ ) and ( [4.16| ), by taking the complex conjugate for the Yukawa couplings and 
replacing Pr with P^. It may not be very convenient to present analytic expressions for 
the CP- violating observables 6n,, defined in Eq. ( [4.10|) , in a rather compact form, although 
their derivation from Eqs. ( |4.19| ) and ( |4.20| ) is quite straightforward. Instead, we shall 
compare the numerical results obtained from our resummation approach with those found 
with different methods. 

Before we proceed with our numerical analysis, it is crucial to take the out-of equi- 
librium constraints on heavy neutrino decays into account. Detailed study of the latter 
will be performed by solving numerically the Boltzmann equations in Section 6. To a good 
approximation however, Sakharov's third condition imposes a lower bound on the lifetime 
of the decaying heavy Majorana neutrino, l/T^., which may qualitatively be given by the 
inequality 

r^,(T = m^J ~ 2H{T = mNj, (4.21) 
where H{T) is the Hubble parameter 



H{T) = 1.73 (^y^ 



M, 



Planck 



(4.22) 



In Eq. ( 4.22 ), g^, ^ 100 — 400 represents the number of active degrees of freedom in usual 
extensions of the SM and Mpianck = 1-2 10^^ GeV is the Planck mass scale. Considering 
the total decay width of the heavy neutrino, = \hii\'^mNj{87i), the out-of equilibrium 
constraint in Eq. ( [4.21| ) yields 



\hiA^ ~ 7.2 X 10-^^ 



1 TeV 



(4.23) 



In order to have the baryon-to-entropy density ratio in the observed ball park, i.e., Yb ~ 
10^^'^ ~ \6i\f-\/g^: , one must allow for CP asymmetries I^atJ of order 10"^ — 10~^. This is 
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practically independent of the heavy neutrino mass, for masses mTv, ~ 1 TeV, provided the 
out-of equilibrium constraint on the Yukawa coupling in Eq. ( |4.23| ) is fulfilled. 




Fig. 4: Numerical estimates of CP asymmetries in scenario I. 
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Fig. 5: Numerical estimates of CP asymmetries versus mN^/iTiNi — 1 in scenario II. 



We shall give numerical estimates of CP asymmetries in two heavy-neutrino scenarios, 
which are in compliance with the out-of-equilibrium limits derived above. For our illustra- 
tions, we analyze models, in which the two-right handed neutrinos mix actively with one 
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lepton family I only. Despite their simplicity, such models exhibit all the essential features 
of CP violation through heavy-neutrino mixing. Specifically, we consider the following two 
scenarios: 

1. mjvi = 10 TeV , hn = 10"^ hn = 10-^(1 + i) , 
II. m^vi = 10^ TeV, hn = 10"^ ha = 10-^(1 + i) , (4.24) 

and assume that N2 is always heavier than A^i, i.e., itlni ^Tn^,^. 

In Fig. 4, we display numerical estimates of the CP asymmetries defined above as 
a function of the parameter = ni]\i.^/mN-^ — 1 for the scenario I. We have divided the 
range of values for the parameter into two regions: The first region is plotted in Fig. 
4(a) and pertains to the kinematic domain, where resonant CP violation due to heavy- 
neutrino mixing occurs. The second one, shown in Fig. 4(b), represents the kinematic 
range, far away from the resonant CP-violating phenomenon. The dotted line in Fig. 




Fig. 6: Numerical estimates of CP asymmetries as a function of mN^/niN^ — 1 in scenario II. 
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4(a) gives the prediction of En^, obtained from Eq. ( ^.12| ) in the conventional perturbation 
theory. Obviously, e^^^* diverges for sufficiently small values of xtv, e.g., xn < 10^^^. If 
resummation of the relevant fermionic self-energy graphs is considered, the prediction for 
Eni is given by the dashed lines in Fig. 4, which shows a maximum for xn ~ 10~^^. In 
such a case, CP violation may resonantly increase up to order of unity [0,^ . As has been 
mentioned above, the physical regulating parameter of the singularity in e^^^* is the finite 
width of the heavy neutrino N2, which arises naturally within our field-theoretic approach. 
Thus, the condition for resonant enhancement of CP violation reads: 

rriNi - '^Ni ~ ±^22"^Af2 = and/or AumN^ = ■ (4-25) 

Clearly, Fig. 4(a) satisfies the above condition. However, the magnitude of the CP asym- 
metries is governed by the expression 

which is always 6c p < 1- Evidently, both scenarios I and II given above represent maximal 
cases of CP violation with 6cp = 1- Therefore, results for any other model may readily be 
read off from Figs. 4, 5 and 6, by multiplying them with the appropriate model-dependent 
factor 6c p. The solid line in Fig. 4 gives the numerical estimate for the CP- violating param- 
eter 6n in Eq. ( [4.11|) , where e'-type contributions are included. The latter are very small in 
this scenario, so as to potentially account for the BAU, e.g., ~ 10^^^. Furthermore, it 
may be important to stress that 6^ vanishes in the CP-invariant limit xat ^ 0, as it should 
be on account of Eq. ( p.8|) . 

In Figs. 5 and 6, we give numerical estimates of the CP asymmetries in the scenario 
II. The difference of this model with the scenario I is that the e'-type effects may not be 
negligible in the off- resonant region, as can be seen from Figs. 5(a) and 6. In particular, 
for values of the parameter x^ < 10^^^ or x^v > 1, the individual e'^^- and E'j^^-tjpe 
contributions may dominate over those of the e type. Models with x^ > 1 have extensively 



been discussed in the literature |]TT],|12|,|T3|,|T4[ . Numerical estimates for such models are 
displayed in Fig. 6. Our attention will now be focused on the domain with xn < 10^^. In 
Fig. 5(a), we observe that e'j^^ and e'j^^, represented by the dotted lines, do not vanish in the 
CP-invariant limit xn — > 0, as opposed to Eni- As a consequence, the CP asymmetry ^tvi 
in Eq. (|4.1CI| ), in which both ENr and E'j^^-tjpe terms are considered within our formahsm, 
does not vanish either. The reason is that the physical CP-violating parameter in this 
highly degenerate mass regime for A^i and N2 is the observable 6n defined in Eq. ( |4.11|) . In 
fact, the quantity 6n shares the common feature with e^-^ and tends consistently to zero 
as xat ^ 0. This fact must be considered to be one of the successes of our resummation 
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approach. Again, CP violation is resonantly amplified, when the condition in Eq. ( ^4.25| ) 
is satisfied, as can be seen from Fig. 5(b). Finally, we must remark that —6n has a zero 
point and eventually becomes negative for xn ^ 1, as is plotted in Fig. 6. Nevertheless, 
this result should be viewed with great caution. The actual reason is that the effect of the 
different dissipative Boltzmann factors multiplying the decay rates of the heavy Majorana 
neutrinos A^^i and N2 must be considered in the definition of 6n- in Eq. ( [4.111) . These 
phenomena will be taken into account in Section 6, in which we shall write down and solve 
numerically the relevant Boltzmann equations. 

Our computation of the CP asymmetries has been carried out in the physical basis, in 
which the heavy neutrino mass matrix is non- negative and diagonal. One might, however, 
raise the question whether our analytic results would have been modified if we had chosen 
a different basis other than the mass basis. The best way to address this question is to 
discuss first how our expressions would change under a basis transformation. Let us assume 
a model with two-right handed neutrinos for simplicity and imagine that we wish to perform 
our calculations in a basis, in which the two heavy Majorana neutrinos, rii and ^2, say, span 
a non-diagonal mass matrix. The heavy Majorana neutrinos ni and ^2 are then related to 
the neutrinos A^^i and A^2, defined in the mass basis, through the unitary transformation: 

The 2x2 unitary matrix relates the physical and non-diagonal mass matrices in the 
following way: 

M" = ?7"^M"?7". (4.28) 

It is now helpful to see how the various kinematic parameters transform under the action 
of [/". In particular, the inverse propagator matrix S'^^(^) undergoes a change depending 
on S:^^{;^), which is given by 

Sn\^) = {U^'^Pr + U^^Pl) S-\^) {U'^Pr + t/"*PL) . (4.29) 
This in turn implies that 

Sn{]^ = {U'^^Pr + U'^^Pl) SM {U'^^Pr + U'^Pl) , (4.30) 

where Sn{j^) and S'„(/i) are the 2x2 propagator matrices, evaluated in the two different 
heavy neutrino bases. Finally, the Yukawa couplings h = {hn, hi2), defined in the mass 
basis, are related to the Yukawa couplings k"- = {hfi, K12) of the non-diagonal basis, via 
the unitary rotation: h = h^U^. 
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Given the afore-mentioned transformation of the Yukawa couphngs and that of the 



resummed propagators in Eq. ( 4.30 ) under a basis rotation, it is not difficuh to show that 



the propagator part of the S-matrix amphtude for the process, e.g., I ^ ^~^X 

^ l^X~) oc TT[hPRS^{^)PRh^] , 
is rotational invariant under U^, i.e., 

TT[hPnS^i^)Pnh^] = Tilh^PnSMPRh''^' 



1 



where the trace should be taken over the product of spinor and flavour matrices as well. In 
fact, this avenue has been followed in If decays of particles are considered however, such 
as heavy Majorana neutrino decays, the above flavour-basis invariance is not manifest. The 
reason is that, within the LSZ formalism, the S-matrix amplitude is obtained by amputating 
the external legs of the Green functions with inverse propagators defined in the mass basis 
and any basis transformation will affect the S-matrix expression for the decaying particle. 

Therefore, it is essential to introduce a method that always makes reference to the 
decay of the neutrinos in the diagonal mass basis. To accomplish this, one has to truncate 
the Green function, e.g., T^^'^Sn^'f), for the decay rii — ^ /'x^; i^i the following scheme: 

S,uM = Tr''SnmU''*PR + U''PL){Zl''^^ Pr + Zr^''^ PlUSm)-^{^)uM . (4.31) 

where = Z*^ due to the Majorana property of the heavy neutrinos and summation over 
not displayed indices is implied. Note that T^'^'^ in Eq. ([4.31|) is calculated in the non- 
diagonal basis. We must remark that our field-theoretic approach to solving the problem of 
fiavour-basis invariance leads to non- vanishing CP asymmetries proportional to the fiavour- 
basis independent combination of CP invariance given in Eq. ( p.8|) . In particular, in the 
limit (mjva — ""^ati) — > 0, the vanishing of in Figs. 4(a) and 5(a) explicitly demonstrates 
that our numerical estimates describe genuine effects of CP violation. 



5 Low-energy constraints 



If heavy Majorana neutrinos are not much heavier than few TeV these novel particles 
may then be produced at high-energy ee [^, ep |^, andpp colliders whose subsequent 
decays can give rise to distinct like-sign dilepton signals. If heavy Majorana neutrinos are 
not directly accessed at high-energy machines, they may have significant non- decoupling 
quantum effects on lepton-flavour- violating decays of the Z boson [|^,^, the Higgs particle 
{H) 1^ , and the r and /i leptons . Their presence may cause breaking of universality in 
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leptonic diagonal Z-boson |^8| and tt decays |^ or influence the size of the electroweak 
oblique parameters S, T and ?7 |50]. In fact, there exist many observables |51], including the 
r-polarization asymmetries, neutrino-counting experiments at the CERN Large Electron 
Positron Collider (LEPl) or at the Stanford Linear Accelerator (SLC), to which Majorana 
neutrinos may have sizeable contributions. However, if the out-of-equilibrium constraints 
are imposed on all lepton-families of the model (c/. Eq. ( [4.23|) ), all these non-SM effects 
mentioned above are estimated to be extremely small at the tree or one-loop level. For 
example, typical tree-level terms breaking charged-current universality in vr decays due to 
heavy Majorana masses depend linearly on the ratio p/m^. and are therefore very 
suppressed. Moreover, one-loop induced flavour-changing decays of the Z boson are also 
very small of order 

\hiiv\'^ i~ 



4vr 



Ml 



< 



10 



-15 



X 



1 TeVx2 

niN,. ^ 



(5.1) 



at the amplitude level, with being the SU(2)j;^ electroweak fine structure constant. 

It is known that below the critical temperature Tc of the first-order electroweak phase 
transition, the Higgs boson acquires a non- vanishing VEV ^ . This leads to non-zero light 
neutrino masses, which may be obtained from the five dimensional operator 



i=l 



mN,, 



(5.2) 



However, for the scenarios I and II, the light-neutrino masses generated at temperatures 
T <^T(. are much below the cosmological limit, i.e., rrii, <^ 10 eV. Obviously, all the afore- 
mentioned limits cannot jeopardize the viability of our minimal new-physics scenarios. 
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Fig. 7: Typical two-loop diagram contributing to the EDM of electron. 

New-physics interactions may also give large contributions to the EDM of electron, 
whose experimental upper bound is {de/e) < 10~^^ cm |]52|. In particular, this bound is 



29 



crucial, since it may impose a constraint on CP-violating operators similar to those that 
induce CP violation in heavy Majorana neutrino decays. In our SM extension with right- 
handed neutrinos, the contribution to the EDM of electron arises at two loops [Q. A 
typical diagram that gives rise to an EDM for the electron, 

U = ^e(^)ea^,75eaM^ (5.3) 



is shown in Fig. 7. Following the semi-quantitative prescription in ||53[ by making naive 
loop counting for rriN^, mjq^^ Mi^, we find 

V le 2e K 21) ^g^2 (^^, + ^^J' ^ Mw ^ 

in-24 ^ c- fu u* u u*\ ^N^mN^jmlj - m% ) /m^,x 
~ 10 cm X Qmih.h^^hM 2 I J,2 ^'^ix^J' ^^-^^ 

where the remaining factor depending on the masses of the heavy Majorana neutrinos only 
in Eq. ( p.4|) is usually less than one. It is then obvious that the above EDM limit may be 
important, if \hii\ > 0.1, i.e., for super-heavy Majorana neutrinos with mjv. > 10^^ TeV. 
On the other hand, stability of the Higgs potential under radiative corrections requires 
that \hii\ = 0{1) [0. Nevertheless, the EDM bound derived above gets less restrictive 
when the mass difference between A^^i and N2 is much smaller than the sum of their masses, 
i.e., {rriM^ — m]\f^)/{mNi +^N2) ^ 1; and is practically absent for values of the parameter 
xjv < 10"'^. As a consequence, in our analysis, we have only considered scenarios with 
l^iil; 1^/2! ~ 10~^, on which the EDM constraint in Eq. ( ^.4|) is still weak. 

In the above discussion of low-energy effects, we have assumed that the out-of equi- 
librium constraint given in Eq. ( [4.23| ) applies to all lepton families. As we have seen in 
Section 4, one lepton family is sufficient to get the CP asymmetries required for the BAU. 
Furthermore, sphaleron interactions preserve the individual quantum numbers B/3 — L^, 
e.g., B/3 — L^.. Also, an excess in Lg will be converted into the observed asymmetry in 
B. Most interestingly, the so-generated BAU will not be washed out, even if operators 
that violate and Lr are in thermal equilibrium provided that possible Le — L^- and 
Le — LT--violating interactions are absent For instance, this can be realized if there are 
two right-handed neutrinos that mix with the electron family only and produce the BAU, 
and the remaining isosinglet neutrinos strongly mix with the /x and r families. This class 
of heavy Majorana neutrino models may predict sizeable new-physics phenomena, which 
have been mentioned above and can be probed in laboratory experiments. 
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6 Boltzmann equations 



In this section, we write down the relevant Boltzmann equations Il3|j^ j55|j56[| , which deter- 



mine the time evolution of the lepton-number asymmetries. Then, we solve these equations 
numerically and present results for the expected BAU within the two different democratic- 
type scenarios I and II discussed in Section 4. Lastly, we give estimates of the finite tem- 
perature effects, which may have an impact on the resonant phenomenon of CP violation 
due to mixing. 

The Boltzmann equations describing the time evolution of the lepton asymmetry for 
a system with two heavy Majorana neutrinos are given by |]55| , [T2|1 

%^^ + 3ffn„. = (6.1) 



—- + mriL = [Sm, - ij _ ^ 7iv, - -ig7., (6.2) 

where 71^., rii = ni — nj are the densities of the number of Ni and the lepton-number 
asymmetry, respectively, and n^. and nf*^ are their values in thermal equilibrium. The 
Hubble parameter H = (dR/dt)/ R determines the expansion rate of the Universe and also 
depends on the temperature T, through the relation in Eq. (|4.22| ). In Eqs. ( |6.1| ) and ( |6.2|) , 
•jNi and 7o- are the collision terms given by 

Kiim%jT) 
K2{m%jT) 

T 



7a 



fOO 

/ dss^/^Ki{^/T)a'{s). (6.4) 
Jo 



Here, Ki{z) and K2{z) are the modified Bessel functions defined in In addition, 
and cr'{s) are respectively the usual T = expressions for the total decay- width of Ni and 
the cross section of the 2^2 scatterings, involving the L and $ states, which are taken 
here to be massless. The latter comprises the scatterings, i.e., L^"^ and its CP- 

conjugate process L^^ L"^^. In fact, the cross section a'{s) is calculated by subtracting 
all those real intermediate contributions that have already been taken into account in the 
direct and inverse decays of heavy Majorana neutrinos IQ. Therefore, 7o- may be regarded 
as an additional CP-conserving depletion term, which can be shown to be of order /ifj in the 
Yukawa couplings, i.e., one formally finds that 7o- ~ 7^. in the narrow width approximation 



m 



In writing down Eqs. ( |6.1| ) and (|6.2|) , several applicable assumptions have been 
made, which are also reviewed in Ref. ^5|. First, we have considered the Friedmann- 
Robertson- Walker model in the non-relativistic limit. Second, we have adopted the 
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Maxwell-Boltzmann statistics, which is a good approximation in the absence of effects 
that originate from Bose condensates or arise due to degeneracy of many Fermi degrees of 
freedom. Third, we have assumed that the lepton and Higgs weak isodoublets, L and $, 
are practically in thermal equilibrium, and neglected high orders in n^./nf'^ and 5^1- In this 
context, it has also been assumed that the different particle species are in kinetic equilib- 
rium, i.e., the particles may rapidly change their kinetic energy through elastic scatterings 
but the processes responsible for a change of the number of particles are out of equilibrium. 
These out-of equilibrium reactions are described by the Boltzmann equations (|6.1|) and 
( |6.2|) . Finally, there may exist additional contributions to the Boltzmann equations fl^ , 
coming from processes, such as NiL — > $* ^ QitR, NiQi Ltji, where Qi {i = 1,2,3) 
denotes the usual quark isodoublets in the SM. These reactions as well as those of the 
kind $<l'^ —>■ LL^ are still very weak to wash out the BAU generated by the direct heavy 
Majorana neutrino decays as long as the out-of equilibrium constraint on the Yukawa 
couplings in Eq. ([4.23|) is imposed. Hence, we have neglected these small depletion terms. 

Before we evaluate numerically the Boltzmann equations written above, it will prove 
helpful to make the following substitutions: 

^ ~ T ' 2HiT) 2H{x = 1) ' ^ ^ 

which is a good approximation for the radiation dominated phase of the Universe. We 
assume the heavy neutrino mass hierarchy rriNi ^ ttin^ for the two- right handed neutrino 
scenarios I and II, given in Section 4. Furthermore, we introduce the quantities Y]\f. = rij^Js 
and Yl = ul/s, where s is the entropy density. In an isentropically expanded Universe, 
the entropy density has the time dependence s{t) = const, x R~^[t) and may be related to 
the number density of photons, n^, as s = g^^n^, where is given after Eq. (|4.22| ). For our 
discussion, it will be more convenient to define the parameters 

^ ^ Ki{x) Tn, ^ K2{x)Ki{^x) Fjv, 

K2{x) H{x = 1) ' ^ K,{x)K2{ix) F^, ' ^ ' ^ 

with ^ = mATj/mTVi- Making use of the above definitions and relations among the parame- 
ters, we obtain the Boltzmann equations for the new quantities Iat^, Yatj and 1^, viz. 



dx 



{Y^,-Y^l)Kx\ (6.7) 



dY, 



- -{Y^,-Y^l)^Kx\ (6. 

1 



dx 
dYL 



dx 



g.YiX^Nl— Kx\ (6.9) 
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In our numerical analysis, we shall neglect the Yukawa coupling suppressed term in Eq. 



Y^lix) = — I dz z y/z^ - x2 e-" = —x^K2{x), (6.10) 



(3.9), which is proportional to 70-, since 70- <^ Moreover, the heavy-neutrino number- 
to-entropy densities in equilibrium Y^^(x) are given by 0] 

3 roo 3 

— / dz z ^/ z'^ — x'^ e"^ = — 

Qsf Jx 8(7* 

and Y^^{x) = F^^(^x). The differential equations (|6.71 )- (|6.9| ) are solved numerically, using 
the initial conditions: 

YnM = YnM = Y^liO) = Y^liO) and rL(0) = 0. (6.11) 

These initial conditions merely reffect the fact that our Universe starts evolving from a 
lepton symmetric state, in which the heavy Majorana neutrinos are originally in thermal 
equilibrium. After the evolution of the Universe until temperatures much below mjvj, a 
net lepton asymmetry has been created. This lepton asymmetry will then be converted 
into the BAU via the sphalerons. During a first order electroweak phase transition, the 
produced excess in L will lead to an excess in B, which is given by 

8Ng + ANh ^ 1^ 
= 22N, + 13Ar^ "^"^ " " 3 ' ^''''^ 

where Yb = ub/s, Ng is the number of generations and A^^^ is the number of Higgs doublets. 
The observed BAU is Y^^^ = (0.6 — 1) x 10^^° [^, which corresponds to an excess of leptons 



_y^bs 10^9 — lO^^''. In the latter estimate, other alternatives for generating the BAU are 
also considered, which may arise from the conversion of an individual lepton asymmetry 
5^ only, e.g., Le, into the BAU. 

In Fig. 8, the observed range for Yl, Y^''*, is indicated with two confining horizontal 
dotted lines. Furthermore, we display our numerical estimates of Yl( OC^ clS db function of 
the parameter x = itlni/T, for selected heavy Majorana neutrino scenarios, stated in Eq. 
( [4.24| ). Specifically, Fig. 8(a) shows explicitly the dependence of Y^ on x, for the three 



different values of the parameter a; at = niN^/mNi — 1 = 10^^, 10^^ and 10~^° in scenario I. 
In this scenario, the lightest heavy Majorana neutrinos A^i has a mass rriN^ = 10 TeV and 
the values of the Yukawa couplings are hn = 10^^ and hi2 = 10^^(1 + i). The parameter 
xat is a measure of the degree of mass degeneracy for Ni and N2. For comparison, it 
is worth mentioning that the degree of mass degeneracy between Kl and Ks is of order 
10^^^, which is by far smaller than the one considered here. Since e'-type CP violation is 
very small in scenario I, as has already been discussed in Section 4, one has to rely on CP 
violation through heavy Majorana neutrino mixing. We find that for small heavy neutrino 
mass splittings determined by x^, with xn being in the range between 10~^ and 10~^, a 
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Fig. 8: Lepton asymmetries for selected heavy Majorana neutrino scenarios. 



sufficiently large lepton (baryon) asymmetry can be generated. The significance of our e- 
type CP-violating mechanism may be seen from the fact that in democratic-type scenarios, 
i.e., in models with all Yukawa couplings being of the same order as those considered here, 
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heavy Majorana neutrinos with masses as low as 1 TeV can still be responsible for the 
excess of baryons, found by observational measurements. Of course, for larger values, 
e.g., xn > 10^®, the BAU is getting much smaller than Y^^*. Furthermore, we have also 
checked that CP violation and hence BAU vanishes in the limit xn — > 0, i.e., when the 
two OS renormalized heavy-neutrino masses Ni and N2 are exactly equal, as it should be 
on account of the CP invariance condition in Eq. (|2.8| ). 

Fig. 8(b) gives numerical estimates of Yl as a function of x, for the scenario II. In 
this model, we have chosen niNi = TeV, and hn = 10"^ and hi2 = 10^^(1 + i). We 
also present results for three different values of the parameter xjy, x^ = 0.1, 1 and 10. 
In this large-mAT^ scenario, a high degree of degeneracy for A^^i and N2 is not required in 
order to get sufficient CP violation for the BAU. In fact, the e- and e'-type contributions 
to the decays of heavy Majorana neutrinos are of comparable order and should both be 
taken into account. Again, one finds numerically an appreciable excess of leptons, Yl, 
within the observed range Y^^'^. For the scenario with x^ = 10, we obtain positive values 
for Yl up to x ~ 0.5. This small excess of leptons is rapidly erased by the Ni heavy 
neutrino decays, which are almost in thermal equilibrium. At lower temperatures, i.e., 
for X ^ 0.5, the heaviest heavy Majorana neutrino N2 gets decoupled from the system, 
and out-of equilibrium A^^i decays will eventually produce a non-zero value for Y^ at the 
observable level. Since the e- and e'-type contributions are formally of order hf- in this 
highly non-degenerate scenario, other collision terms of order hf^ may also be significant, 
such as the scatterings ^L*-^ and LL'~^ |112|,p8|. However, the inclusion of 



these additional effects will not quantitatively affect our numerical results much, as long 
as the constraint in Eq. (^4.23|) is valid. 

In our numerical analysis presented above, we have not taken into account other ef- 
fects, which might, to some extend, affect the resonant condition, given in Eq. ([4.25|) . Apart 
from the intrinsic width of a particle resonance, there may be an additional broadening at 
high temperatures, due to collisions among particles. Such effects will contribute terms of 
order /i| to the total Ni widths and are small in general P,p9|. Of most importance are, 
however, finite temperature effects on the T = masses of the particles. Since SM gauge 
interactions are in kinetic equilibrium in the heat bath, they can give rise to thermal masses 
to the leptons and the Higgs fields ||60| , |6T| , |62| . These thermal masses are given by 

^ = ^ (3/ + 9") - 0.044 
M|(T) ^ T2 



T2 2d(l-^), (6.13) 

where g and g' are the SU(2)i and U(l)y gauge couplings at the running scale Mz, respec- 
tively, and d = [2M^ + M| + 2mf + M^]/{8v^). The critical temperature calculated at 
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one loop may be obtained from |62 
1 



Ml-^, (2M^ + Ml - Amt) - (2M^ + M^f] . (6.14) 



Although the isosinglet heavy neutrinos do not have tree-level couplings to the SM gauge 
bosons, the difference between the thermal mass of Ni and its respective zero-temperature 
mass will proceed through Yukawa interactions |60], 



I.e. 



m 



,(T) - ml,SO) _ 1 



"'"^2 • = 1^1""!' • (6.15) 

Such a T-dependent shift in the masses of is very small and may be safely neglected in 
the mass difference (T) —m^^ (T), which enters the analytic expressions for resonant CP 
violation (see, e.g.., Eqs. (|4.12|) and (f4.13|) ). Making now use of Eqs. (|6.13|) and (|6.14|) , the 



authors in fe^ find that 0.5 < Af$(T)/T < 2, when Mh varies from 60 GeV up to 1 TeV, 
for T :^T^^ 200 GeV. In particular, one obtains M^{T)/T ~ 0.6, for Mh < 200 GeV. 
In this Higgs-mass range, the effective decay widths of the heavy neutrinos rAr^(T) will be 
reduced relative to rAr^(O) by 70% - 80%, because of considerable phase-space corrections 
|6T| . If Mh > 350 GeV, M^(T)/T is getting bigger than one, which signals the onset of a 



non-perturbative regime and pure perturbative methods may not be sufficient to deal with 
very heavy Higgs bosons. For temperatures T near to the critical temperature Tc, M$(T)/T 
will be smaller because of the suppression factor (1 — T^/T'^) in Eq. (|6.13| ). It appears that 



low-scale leptogenesis is less affected by finite temperature effects, even though one can 
always choose larger Yukawa couplings to enhance rfyf.{T) in the light-Higgs scenario. In 
either case, the resonant phenomenon of mixing-induced CP violation plays a crucial role 
to generate sufficiently large CP asymmetries. 



7 Conclusions 

We have studied the impact of the e- and e'-tjpe mechanisms for CP violation on generating 
the excess of baryons detected in the Universe. As for the scenario of baryogenesis, we have 
considered that out-of-equilibrium L- violating decays of heavy Majorana neutrinos produce 
an excess in L, which is converted into the observed asymmetry in B, through the B + L- 
violating sphaleron interactions. In Section 2, we have described minimal extensions of 
the SM with right-handed neutrinos, which can predict nearly degenerate heavy Majorana 
neutrinos without resorting to a fine tuning of the mass parameters. Such models may 



naturally occur in certain subgroups of SO(IO) pSf or Eg theories |2^. In a one-lepton 



family model, the presence of two right-handed neutrinos is sufficient to give rise to the 



36 



non-trivial CP-violating combination of Yukawa couplings, given in Eq. ( p.SI ). As has been 
demonstrated in Section 4, our physical CP asymmetries depend indeed on this CP-odd 
invariant. In addition, particular emphasis has been laid on the renormalization of the 
Yukawa couplings in the minimally extended model. 

In the conventional perturbation theory, the wave-function amplitude becomes singu- 
lar, whenever the degenerate limit of the two mixed heavy neutrinos is considered. Several 
effective methods have been proposed to solve this problem, such as diagonalizing the ef- 
fective Hamiltonian of the two-heavy-neutrino system |1T6| , ^ . The results obtained with 
these methods show a resonant enhancement of CP violation, when the two heavy neu- 
trino masses are getting closer. Such a resonant CP-violating phenomenon is in line with 
earlier studies in |p2|. Unfortunately, the methods based on diagonalizing the effective 



Hamiltonian are not analytic, if the effective Hamiltonian itself is not diagonalizable p5 



Therefore, it is important to compare the results obtained for the CP asymmetries with 
a more rigorous field-theoretic approach. In Section 3, we have extended, in an effective 
manner, the resummation formalism for particle mixing, which was applied to scatterings 



in to that of the decays of particles. The resummed decay amplitudes possess all the 
desirable field-theoretic properties and exhibit an analytic behaviour in the mass degenerate 
limit. 

In Section 4, we have used the afore-mentioned field-theoretic approach to perform a 
systematic analysis of the e and e' types of CP violation in the L-violating decays of heavy 
Majorana neutrinos. For illustration, we have considered minimal extensions of the SM with 
two right-handed neutrinos. We have found that e-type CP violation is resonantly enhanced 
up to order of unity, if the mass splitting of the heavy Majorana neutrinos is comparable to 
their widths, as is stated in Eq. ( [4. 251) , and if the parameter 5c p defined in Eq. (|4.26|) has 
a value close to one. In our view, these two necessary and sufficient conditions for resonant 
CP violation of order unity constitute a novel aspect in the leptogenesis scenario, which 
have not been pointed out in their most explicit form before. Taking full advantage of the 
mechanism of resonant CP violation, one may consider scenarios with nearly degenerate 
heavy Majorana neutrinos at the TeV scale and all Yukawa couplings being of the same 
order, which can still be responsible for the BAU. In this kinematic range, the 5'-type 
contributions are extremely suppressed. Of course, the higher the isosinglet Majorana mass 
is, the less the above degeneracy is required in order to get sufficiently large CP violation. 
However, even in the weak- mixing limit, i.e., tun^ <^ t^N2i ^-type CP violation is equally 
important with the e'-type one and therefore should not be ignored |T5|. In contrast, only 
e-type CP violation is important in the strong mixing regime, i.e., when rriN-^^ —771^2 ^'^n,- 
Another alternative of having sufficiently large CP violation for TeV leptogenesis, with 
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mjv^ = 0{1) TeV, is to assume an hierarchic pattern for the heavy Majorana masses and 
the Yukawa couphngs, e.g., m^.^ ^ "^a^i hi2 ^ hn. Such scenarios were thoroughly- 



investigated in ||T^,r^|, and therefore we have not repeated this analysis here. On the 
other hand, a wide range of heavy neutrino masses, Tc ~ mjVi ~ 10^ TeV, is still able to 
account for the BAU, even if all Yukawa couplings are of the same order. We can hence 
conclude that the two CP-violating mechanisms under consideration are, to a great extend, 
determined from the flavour structure of the neutrino mass matrix and the flavour hierarchy 
of the Yukawa couplings in the model. 

In Section 5, we have presented estimates of low-energy constraints on our minimal 
model, coming mainly from the electron EDM. The two-loop EDM bound derived is found 
to be not very severe in order to rule out the leptogenesis scenario and is practically absent 
if the heavy Majorana neutrinos are nearly degenerate. In Section 6, we have briefly 
discussed the implications of finite temperature effects for our resonantly enhanced CP 
violation. Although temperature effects on the heavy Majorana neutrino masses are very 
small, because isosinglet neutrinos interact quite feebly with the Higgs fields in the thermal 
bath, the final leptons and Higgs fields acquire appreciable non-zero thermal masses. As 
a consequence, there will be a reduction of the widths of the heavy Majorana neutrinos, 
relative to their respective values at T = 0, due to a phase-space suppression. The size 
of the reduction of the Ni widths depends crucially on the zero-temperature Higgs-boson 
mass Mh- For Mh < 200 GeV, the effective decay width of N^, rjv-(T), is estimated to be 
smaller than r7v.(0) by 80% at most. This will roughly lead to an 80% decrease of the CP 
asymmetries, calculated at zero temperature. In this respect, the resonant phenomenon of 
CP violation through mixing of heavy neutrinos plays a very important role for leptogenesis, 
since e-type CP violation can still be large for heavy-neutrino mass differences comparable 
to the effective decay widths T^XT). Finally, further support for the viability of the 
resonantly enhanced CP-violating phenomenon is obtained from solving numerically the 
Boltzmann equations. Numerical estimates reveal that Eg-type scenarios, which naturally 
predict a certain degree of degeneracy between the heavy Majorana neutrinos, are able 
to account for the present excess of baryons. In conclusion, even if all heavy Majorana 
neutrinos have masses as low as 1 TeV and all couple to the lepton and Higgs isodoublets 
with a universal Yukawa strength, they can still be responsible for the BAU observed in 
our epoch, by means of the resonant mechanism of CP violation presented in this paper. 
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A One-loop analytic expressions 



In this appendix, we list the analytic expressions for the one-loop self-energies of the 
Higgs and fermion fields as well as the one-loop vertex couplings x^^^ii X^^i^i ^i-nd HuiNi. 
Detailed discussion of mixing renormalization for Dirac and Majorana fermion theories may 



be found in and will not be repeated here. Instead, we present the relations between 
the wave-function CT's and unrenormalized self-energies. Our analytic results will be 
expressed in terms of standard loop integrals presented in adopting the signature for 
the Minkowskian metric g^^i, = diag(l, —1, —1, —1) (see also Appendix A of Ref. [0]). 

The Feynman rules used in our calculations may be read off from the Lagrangian ( p.4| ). 
We first compute the Higgs self-energies X^X^ ^ind HH, shown in Fig. l(d)-(f). They 

all found to be equal, viz. 



n 



X X 



riL riR 

EE 

1=1 i=l 



\hii 



(A.l) 



Thus, the universality of the divergent parts of the wave functions SZ^- , SZ^o and SZfj is 
evident, if one calculates SZ^^"" = — 9f?en'^™(0) from Eq. ( |A.1| ), for all field components of 
the Higgs doublet $. 

The individual contributions to the one-loop fermionic transitions, /' I, z/// ui 
and Nj Ni are displayed in Figs. 1(g), 1(h) and l(j), respectively. Explicit calculation 
of the fermion self-energy transitions gives 



E^^^.5l(p^<,o) 

1=1 



-E 



1 



(A.2) 



^ 327r2 



+ mNAKhi'iPR + hiihi,,PL)[Bi{p^,m%^,0) - Bi{p^ 
-E ^Pr + huK^ ^Pl){\b,{p\Q,Q) 

+ \b,{p\q,mI) 



(A.3) 



(A.4) 
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Note that the hght-neutrino self-energies contain non- vanishing mass terms in the 

hmit ^ — s> 0. Even though these contributions vanish when 0, they are non-zero at 

T = 0, because Mh 7^ 0, and hght neutrinos may hence receive small radiative masses ||65|| . 
However, the latter are generally controlled by the mass differences of the heavy neutrinos 
and/or the Higgs Yukawa couplings hn. If the range of parameters relevant for the BAU is 
considered, as has been derived in Section 6, possible experimental limits on light neutrino 
masses are estimated to be not very restrictive. 

The wave-function renormalization constants can now be expressed in terms of un- 
renormalized self-energies |38]. Before doing so, we first notice that the one-loop fj —>■ fi 
transitions calculated above between the fermions fi = I, z//, Ni have the generic form 



R 



(A.5) 



where only dispersive parts are considered. If the transitions are between Majorana 
fermions in Eq. ( |A.5| ), one then has the extra properties = and S^(p^) = 

E^(p^). Adapting the results of to our model, we obtain the wave-function CT's 



5Zl 



and, for i ^ j, 



+ 



2 m,- 



rrii 



(A.6) 



5ZL 



mt 



(A.7) 



The wave-function renormalization of charged leptons may be recovered from Eqs. (|A.6|) 
and ( |A.7] ), if one drops all terms depending on S*^ ijp') and its derivative. At this point, it is 
important to remark that there will be additional contributions to our wave-function CT's 
of the Higgs and fermion fields from the gauge sector of the SM . For example, there are 
non- zero thermal mass effects on the particles involved in the loop, which may further break 
the wave-function universality of the different components of the Higgs doublet |^,^. As 
has been discussed in Section 2, this should not pose any problem to the Yukawa coupling 
renormalization, as long as all differences of the kind 5Z^- — 6Zh or SZ^- — SZ^^o are UV 
finite. 

We will now present analytic results for the one-loop corrections to the vertices 
X^PNi, x^^i^i and HuiNi, keeping all M//-dependent mass terms. From Figs. 1(a)- 
(c), we observe that the one-loop coupling x^/^A^j proceeds via a AL = 2 mass insertion 
only, whereas the couplings x^^i^i ^i-nd HviNi can occur through both L-conserving and 
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L-violating interactions. More explicitly, we have 



-iV 



HviN, 



IUiPrUn, ^g^2 ' f^l'ihl'jhlj 



l'=lj=l 
,2 n ™2 



Co(0,0,m^^,0,m^ 0) + Ci2(0,0,m^^,0,m^ 0) 



(A. 



- UiPrUn, { 



V=lj=l 



32^2 



h'l'ihi'jhij 



X 



Co(0,0,m^^,0,m^^.,0) + Co(0, 0, m^^, M^, m^^, 0) 

+Ci2(0,0,m^,,0,m^^,,0) + ^12(0, 0, m^^, M^, m^^ , 0) 
1 

— h„ h:. . 



32\/27r2 



hi'ihlijhij 



m7v^mAr^(Ci2(0,0,m^^,0,m^,,0) - Ci2(0, 0, m^^, M|, m^^, 0)) 
-M^Co(0,0,m^^,Mi,m^^,,0)] }, 



(A.9) 



iUiPRUN, 



hi'ihiijhij 



Co(0,0,m^^,0,m^^,,0) + Co(0, 0, m^^, M|, m^^, 0) 

+Ci2(0,0,m^^,0,m^^,,0) + Ci2(0, 0, m^^, M^, m^^ , 0) 
1 



hi'ih*ii-hij 



X 



miv^rnAr^(Ci2(0,0,m^^,0,m^^,0) - Ci2(0, 0, m^^, M^, m^^, 0) 
M^Co(0,0,m^^,M^,m^^,,0)] }. 



(A.IO) 

As can be readily seen from Eqs. ( |A.9| ) and ( [A.IOD , the part of the couplings x^^i^i ci^nd 
HuiNi that conserves lepton number is UV finite and vanishes identically in the massless 
limit of the Higgs boson. 
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